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ABSTRACT 

We find a significant number of massive and compact galaxies in clusters from the ESO Distant Clusters 
Survey (EDisCS) at 0.4 < z < L They have similar stellar masses, ages, sizes and axial ratios to local z ^ 0.04 
compact galaxies in WINGS clusters, and to z = 1 .4 - 2 massive and passive galaxies found in the general 
field. If non-BCG cluster galaxies of all densities, morphologies and spectral types are considered, the median 
size of EDisCS galaxies is only a factor L18 smaller than in WINGS. We show that for morphologically 
selected samples, the morphological evolution taking place in a significant fraction of galaxies during the last 
Gyrs may introduce an apparent, spurious evolution of size with redshift, which is actually due to intrinsic 
differences in the selected samples. We conclude that the median mass-size relation of cluster galaxies does 
not evolve significantly from z ^ 0.7 to z ^ 0.04. In contrast, the masses and sizes of BCGs and galaxies with 
M^> 4x10"Mq have significantly increased by a factor of 2 and 4, respectively, confirming the results of a 
number of recent works on the subject. Our findings show that progenitor bias effects play an important role in 
the size-growth paradigm of massive and passive galaxies. 
Subject headings: galaxies: clusters: general — galaxies: evolution — galaxies: structure 



1. INTRODUCTION 



High-z studies (as far as z '^ 2.4 ) have found a signifi- 
cant number of massive, passively evolving galaxies (stellar 
mass M* > 10'"Mfr)) with relatively small effect ive radii R^ < 



2kpc ( see, among o thers, iTruiillo et al. 20061 ICimatti et al 
van Dokkum et alj 120081: van der Wei et al. 200? 



2008t 



Saracco et al.ii2009i) . sometimes named superdens^ galaxies 



(SDGs). The general claim by various authors is that lo- 
cal galaxies are three to six times larger in size when com- 
pared to high - z one s, at the same stellar mass. In addition, 
ITruiillo et al.l (120091) found a complete absence of massive, 
old and extr emely compact galaxies in the local universe. 

However, IValentinuzzi et al.l (120101) (hereafter VIO) have 
shown that 22% of local cluster members in the WINGS sam- 
ple with M,^ > 3xlO'"M0 and Esq > SxlO^Mekpc'^ have the 
same characteristics of the high-z SDGs reported in the liter- 
ature by various authors. In the same paper, the authors found 
that selecting galaxies with old stellar populations is equiv- 
alent to selecting the smaller ones, for a given stellar mass. 
Since a large number of galaxies have stopped forming stars 
at relatively low redshift (z < 1 .4), and these tend to be the 
largest, it is not valid to compare high-z passive galaxies with 
all low-z passive ones. To avoid selection effects when mak- 
ing comparisons with passive galaxies at high redshift, one 
needs to select locally those galaxies which were already pas- 
sive at the cosmic time the high-z data correspond to. 

More recently. iTaylor e t al. (2009) revisited the search of 
SDGs in SDSS-DR7 and found a relatively small but signif- 
icant number of SDGs. Following the same criterion used in 
VIO, they find a L3% fraction of SDGs. 

The issue is much debated. iMancini et alj (120101) have an- 

' Regarding physical densities, these galaxies are anyway thought not to 
be extreme (see, e.g.. Bezanson et al. 2009 : Hopkins et al. 2009) 



alyzed a sample of 12 galaxies at 0.5 < z < 1.9 in the Cos- 
mos field, finding masses and sizes compatible with the local 
SDSS ones. Furthermore, by using a set of simulated early- 
type galaxies, they have shown that the low signal-to-noise of 
high-z images can cause measured effecti ve-radii to be lower 
than the intrinsic values. In a recent paper va n Dokkum et alJ 
(|2010) select galaxies with a constant number density at dif- 
ferent cosmic times. They use all galaxies instead of only 
passive ones, and find that galaxies have grown in size by a 
factor of 4 from z '^ 2 to z ~ 0^ 

Even more recently, while ISzomoru et al.l (1201 0*) confirm 
the extr eme compactness of a z = 1.9 galaxy with the HST- 
WFC3, ISaracco et al.l (12010) show that the comoving num- 
ber density of compact ETGs over the volume of about 
4.4xlO^Mpc^ sampled by the GOODS area between 0.9 < 
z < 1.92 is compatible even with the local lower limits given 
in VIO. 

In this Letter we present the results of a search for SDGs in 
the ESO Distant Clusters Survey (EDisCS) at z -- 0.7, and we 
report the comparison of the mass-size relation (MSR) with 
the same relation in WINGS clusters at z '^ 0. We further 
discuss selection effects which may introduce a spurious size 
evolution with redshift if not properly taken into account. 

2. THE DATA 

The high-z cluster sample is extracted from EDisCS, a 
multiwavelength photometric and spectroscopic survey of 
galaxies in 20 fields containing galaxy clusters at 0.4 < 
z < 1 (Whit e et a l. 2005) . We will use a sub-sample of 8 
clusterq^ which have HST-ACS image s for high-precision 
size measurements (iDesai et alJ |2003), and cluster cen- 

^ Their short names found in other EDisCS's papers are CL1138,CL1138a, 
CL1040, CL1216, CL1054-11, CL1054-12, CL1232, CL1354. 
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tral velocity dispersions (adus > 400km s \ < ddus >~ 
700km s~') similar to local WINGS clusters (iHallidav et all 
l2004t iMilvang- Jensen et aTllIOOSt iDesai et al.ll2007l ). Three 
of these clusters have z ~ 0.5, the rest of them have z = 
0.7-0.8. 

G alaxy stellar masses were estimated using the kcorrect 
tool (iBlanton & Roweisl l2007rL that models the available ob- 
served broad band photometry (VRIJK, or BVIJK), fitting 
templates obtained with spectrophotometric models. The stel- 
lar masses are defined as the mass locked into stars, including 
stellar remnants, at any time, using a Kroupa (2001) initial 
mass function (type 2 mass in VIO). Taking into account the 
statistical errors on the mass estimates, the error of the stel- 
lar mass on individual galaxies is of the order ^ 0. Idex, even 
though it has to be taken into account that the scatter (rms) in 
the relation between masses computed with different models 
is typically ^ 0.2dex (for further details, see[ ^itz et al.ll2067l : 
iLonghetti & Saracco..2009:.Vulcani et al1l2010l) 

We use visual morphological classifications from 
IDesai et al.1 (l2007h . 

We measure galax y effective-radii Re with the GIM2D tool 
JSimard et alj|2002l) on the HST images in F814W band, by 
using a single component Sersic fit. The circularized Re is 
determined by numerically integrating the curve of growth of 
the fitted Sersic model, and solving th e equation Flux(< Re) = 
0.5 • Flux(cx)) (for further details, see ISaglia et alj|2010 , sub- 
mitted). The typical random error on the EDisCS's sizes is of 
the order of 20% (Simard et al. 2009). 

We use a mass limited sample of EDisCS spectroscopi- 
cally confirmed cluster members, with stellar masses > 4x 
1O'°M0. This mass limit corresponds to the mass of an ob- 
ject whose observed magnitude is equal to the faint magni- 
tude limit of the spectroscopic survey, with the reddest possi- 
ble color. We correct for spectroscopic incompleteness using 
iMilvang- Jensen et al.l (12008.) completeness functions. 

The local sample examined in this Letter comes from 
the W ide-field Nearby Galaxy-clusters Survey ( iFasano et alj 
120061) . WINGSQ is a multiwavelength photometric and spec- 
troscopic survey designed to provide a robust characterization 
of the properties of galaxies in nearby clusters. 

We use only cluster members of the subset of WINGS clus- 
ters that have an average spectroscopic completeness larger 
than 50% (21 out of 78 clusters), and correct for sp ectroscopic 
incompleteness using the prescriptions given in iCava et al] 
(|2009). These WINGS clusters have redshifts 0.04 <z< 0.07 
and central velocity dispersions 558 < (Tc/us/km s"' < 1368. 

WINGS effective-radii, axial ratios and Sersic indexes 
are measured on th e V-band images with GASPHOT 
dPignatelh et al J 12006 1). an automated tool which performs a 
simultaneous fit of the major and minor axis light growth 
curves with a 2D flattened Sersic-law, convolved by the ap- 
propriate, space-varying PSF. As a measure of galaxy size we 
use the circularized effective radii, calculated in the same way 
it was done for EDisCS sizes (see above). We note here that 
SDG fractions and number densities are updated accordingly 
in this Letter compared to VIO (see next sections). The max- 
imum error on WINGS sizes, ba sed on exte nsive simulation 
runs, is of the order of 10% (see. lPignatelliet al. 2006). 

As a consistency check on sizes, we run GIM2D on one 
representative V-band WINGS cluster image, to compare the 
resulting circularized R,, and Sersic index n of ^ 800 galaxies 



http://cosmo.nyu.edu/mbl44/kcorrect/ 

http: //web. oapd. inaf. it/wings 



with GASPHOT values. We found a systematic difference in 
sizes of 0.033 ± 0.002dex, in the sense that GASPHOT sizes 
are larger than GIM2D ones. This difference becomes larger 
(as far as ~ 0.3dex) for larger galaxies, somehow confirm- 
ing that GIM2D has the tendency to systematically underes- 
timate the sizes of the largest galaxies at all luminosities (see 
ISimard et alj|2002l) . On the other hand, we do not find any 
systematic difference regarding the Sersic index estimate. 

Stellar masses of WINGS galaxies have been determined 
by fitting the optical spectrum (in the range ^ 3600^ ^ 7000 
A), with the spe ctro-photometric model fully described in 
iFritz et all (l2007h . and correcting for color gradients outside 
of the fiber (see VIO). The model derives the integrated spec- 
trum as the combination of stellar populations of 13 different 
ages, allowing dust extinction to vary with the stellar popula- 
tion age and using the single metallicity (either Z=0.05, 0.02 
or 0.004) that gives the lowest y^ fit of the observed spec- 
trum. Although the masses were calculated in two different 
ways we have shown in VIO (and soon in Fritz et al. 2010 
in prep.) that there is no significant systematic offset between 
different methods that could be capable of biasing our results. 

WINGS morphologies are derived from V images using the 
purposely devised tool MORPHOT. We have verified that the 
differences in classification between MORPHOT and an ex- 
perienced human classifier are comparable to the differences 
between two experienced human classifiers (Fasa noet alj 
2010b, in preparation). 

For the sake of comparing the median sizes of the two sur- 
veys, we divide the total sample into four mass intervals, se- 
lected to have a statistically significant number of objects in 
each one of them: 

• BINl: 4x10"' < M,/M0 < 6xl0'° 

• BIN2: 6x10'" < M,/Mq < IxlO" 

• BIN3: 1x10" < M,/Mq < 2x10" 

• BIN4: 2x10" < M./M© < 4x10" 
and will refer to them with the label BIN[l-4]. 

3. EDisCS SUPERDENSE GALAXIES 

In Figure[T]we present the MSR (mass-size relation, bottom 
panel) and the mass-density relation (top panel) of EDisCS 
cluster members with M* > 4xlO'"M0. Colors differenti- 
ate the morphological types (see caption and legend), large 
op en squares are the brightest cluster galaxies (BCGs) listed 
in White et al. 72005). 

In the top panel, the dashed line isolates the EDisCS SDGs 
(larger dots in both panels) with the same density selection 
criteria (E50 > 3xlO^MQkpc"^) used in VIO, above the mass 
completeness limit of this Letter These criteria were cho- 
sen to select galaxies with mass and density ranges similar to 
those of high-z (z > 1 .4) passively evolving galaxies. 

As apparent in Figure [T] we do find a significant number 
of SDGs in the EDisCS sample. Indeed, EDisCS SDGs rep- 
resent 41% of the total cluster population of galaxies more 
massive than M* > 4xlO'"M0. This is an even larger frac- 
tion than VIO found in WINGS local clusters, where 17% are 
SDGs for the mass limits and radii adopted in this Letter A 
decline with time of the SDG fraction in clusters might be 
expected given that a) the "oldest" galaxies in the Universe 
(those who stopped forming stars very early on) inhabit clus- 
ters since very high redshifts, and clusters accrete throughout 
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Figure 1. Circulaiized effective radius Re and mass-density inside Re as a function of stellar mass for all EDisCS spectroscopic members galaxies with M, > 
4x1O"'M0. The different colors mark the morphological type: blue for late-types (later than SOs), green for SOs, red for ellipticals and black for galaxies without 
a classification. Bigger dots highlight the SDGs. Big open magenta squares are the BCGs. The soUd and dashed lines in the bottom panel are the median 
completeness weighted MSRs of EDisCS and WINGS, respectively, obtained excluding the BCGs. Error bars are lower and upper quartiles of the medians. The 
WESfGS mass medians are shifted by O.Oldex in X to avoid overlapping. 



their history galaxies with more extended star formation his- 
tories, and b) as shown in VIO, at any given mass the oldest 
galaxies tend to be the most compact. Therefore the original 
population of old and compact galaxies in clusters get pro- 
gressively diluted by larger galaxies infalling into clusters at 
later times. 

Of the EDisCS SDGs, 41% are ellipticals, 36% are SOs, 
20% are late-type galaxies, and for 4% of them it was not pos- 
sible to assign a reliable visual morphological classification. 
In Table [T] we present the main mean properties of EDisCS 
SDGs. 

In the bottom panel of Figure[T]the median MSR of EDisCS 
galaxies is presented as a black solid line, and compared to the 
dashed WINGS one. As we did in VIO we excluded BCGs 
and galaxies with M» > 4x10"Mq, that are discussed sepa- 
rately below. We find that the median Rg of EDisCS cluster 
galaxies with 4xl0'° < M^/M© < 4x10" is only a factor 
1.18 lower than the WINGS one. Considering separate mass 
bins, the maximum amount of evolution is 1.48 (BIN3, see 
Table |2]i, while in the other mass bins the median sizes turn 



Table 1 

Completeness corrected quantities of EDisCS and WINGS SDGs. EiTors on 
the medians are reported too. 



Quantity 


EDisCS 


WINGS 


SDG fraction 


41% 


17% 


Ellipticals 


41% 


28% 


SOs 


36% 


64% 


Late-type 


20% 


8% 


Unknown morph. 


3% 


- 


Eff. radius (R^) 


1.70 ±0.08 


1.79 ±0.04 


Sersic index (n) 


3.71 ±0.14 


3.21 ±0.09 


Axial ratio {b/a) 


0.59±0.I1 


0.62 ±0.03 


Stellar mass (M*) 


(I.O8±O.O8)x1O"M0 


(l.O2±O.O4)xlO"M0 



out to be in good agreement. 

These numbers may be compared with field studies that are 
includi ng all morphol o gical and spectral galaxy types. Re- 
cently, IWilUams et alj ( 1201 Oh report a size evolution for all 
galaxies from z ~ 0.7 to z ~ 0.04 of a factor of ~ 1.4. In 
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Figure 2. Compaiison of the mass-size relation of EDisCS (top panel) and WINGS (bottom panel). Color coding is the same of Figure[T] The black straight line 
delimits the area above which there are no early-type galaxies in EDisCS. The histogram at the bottom-right represents the luminosity- weighted age distribution 
of WINGS early-type galaxies above (blu dashed line) and below (soHd red line) this line. The black solid and dotted (shifted by -O.Oldex in mass) lines are the 
median MSR for early-type and all types of galaxies in EDisCS, respectively. The dashed black line is the median MSR for WINGS early-type galaxies (shifted 
by +0.()ldex). Error bars are lower and upper errors on the medians. 



Table 2 

Ratios of median WESfGS/EDisCS sizes for the different mass intervals (see 
text). Errors come from the standard error propagation technique. 



WINGS/EDisCS 


BINl 


BIN2 


BIN3 


BIN4 


All galaxies 


1 16+0 23 
'■'"-0.17 


^.2^^^ 


1 4X+«56 


1 nv+o'^ 


Early-type galaxies 


1 7<;+0.19 


^^n^ 


1 f,7+0.67 


1 To+0.95 
^■^'^-0.88 


WINGS Early/EDisCS All 


1 n+o-z 


1.18-- 


1 40+0.57 


1 Ol+O'S 



this regard, we stress that it is hard to draw conclusions from 
this comparison because it is still unclear how the incidence 
of SDGs a nd the evolution of galaxy sizes depend on environ- 
ment (see. lMaltby et al.ll2010l: lRettura et al/'JOlO*). 

So far, we have seen that a considerable fraction of EDisCS 
cluster members are SDGs and that galaxy sizes in EDisCS 
and WINGS, at all mass ranges considered, are rather similar, 
and do not suggest a strong increase in size with redshift. 

The BCGs and the most massive cluster galaxies with M,j> 
4x1O"M0 have to be discussed separately, due to their pe- 



culiar nature and evolution (see, amongst others. Fasano et al] 
l2010al) . Indeed, the EDisCS BCGs have mean mass and size 
of yV^ ^ 4x1O"M0 and Rg ^ 8.5kpc, respectively. In con- 
trast, WINGS BCGs have mean values of M* -- lO'^M© and 
Re ^ 33.6kpc, suggesting that the mean size and mass of 
BCGs have increased by a factor of '^ 4 and ^ 2 between 
z ~ 0.7 and z ^ 0.04, respectively. A lthough this result seems 
in contrast with'Whil ev et alj (|2008), we note that the mass of 
local BCGs in that paper was calculated inside an aperture of 
37kpc, which is approximately the median half luminosity cir- 
cularized radius in V band of our local sample of BCGs. This 
is consistent with a picture where the BCG progenitors in- 
crease their mass via minor mergers in the outer regions , leav- 
ing practically unchanged the dense core (see Hopkin s~et alj 
I2OIOI) . We also note that the size and mass evolution of our 
sample of high-z BCGs with redshift is compatible with the 
observational study of Bernardi (2009) and with the theoreti- 
cal expectations of De Lucia & Blaizot (i2007i) . that predict a 
factor of 3-4 growth in mass between z ~ 1 and z = 0. 

4. SELECTION EFFECTS 
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We have seen that the morphological fractions among 
the EDisCS SDGs are considerably different from WINGS 
SDGs. The latter show a larger fraction of SOs and a cor- 
responding lower fraction of later types. This is expected, 
as many studies have come to the conclusion that a large 
fraction of today's passive early-type galaxies have evolved 
from star forrning late-type galaxy progenitors in clusters 
JDressiereFaDfTW, iFasano et al.ll2000t iPostman et ai]|200l 
iSggianti et al. 2009). 

In VIO we have shown that selecting the oldest cluster 
galaxies means selecting the smallest in size. In the follow- 
ing, we will highlight the biases that can be introduced by se- 
lecting galaxies morphologically, and thus the importance of 
properly taking into account the morphological change too. 
Although the morphological evolution is strictly l inked to th e 
evolution in star formation activity (Poggianti e t alj I2009t) . 
the time scal es can be largely different (Poggianti et al.l ll999t 
iSanchez -Blazquez et al. 2009) and thus become important at 
different cosmic times, in a way difficult to predict. 

In Figure |2] we compare EDisCS (top panel) and WINGS 
(bottom panel) MSRs. Color coding is the same as for Fig- 
ure [1] The black solid and dashed lines show the median 
MSR for early-type galaxies in EDisCS and WINGS, respec- 
tively. Error bars are errors on the medians. The median 
size of WINGS early-type galaxies is a factor 1.53 larger than 
EDisCS's early-types, reaching an offset as large as 1.7 at the 
lowest masses (BINl and BIN2, see Table|2]i. 

At face value, this could be interpreted as an evolution in 
the sizes of individual early-type galaxies. However, we note 
that the largest EDisCS cluster members tend to have late-type 
morphology (some of which are star-forming, ^ 70%, and 
some passively evolving, ^ 30%). We have arbitrarily iden- 
tified a region in the mass-size diagram, above the tilted line 
drawn in both panels, where EDisCS galaxies are large and 
all have late-type morphologies. In contrast, WINGS galaxies 
in this region are mostly (72%) early-types, consistent with a 
morphological evolution. 

A convincing test of this picture is presented in the 
lower right inset of Figure |2] where the distributions of 
the luminosity-weighted ages of WINGS early-type galaxies 
above (blue dashed) and below (red solid) the tilted line is 
shown. The blue-dotted histogram is visibly sharply peaked 
toward lower ages, when compared to the red one. This is 
consistent with a significant fraction of EDisCS large, late- 
type galaxies having turned into large, passive, early-types by 
z-0. 

Let's now focus only on late-types turning in SOs. Practi- 
cally all EDisCS SOs are SDGs (see Figure [U and most of 
them are old (81% have Luminosity-weighted-age > 3 Gyr, 
52% > 5 Gyr); in WINGS, instead, a large number of SOs 
are not SDGs (see Figure |2|i. On the other hand, among 
the WINGS SDGs with SO morphology, only 20% have ages 
lower than the corresponding EDisCS lookback time, i.e. were 
most likely morphologically changed at redshifts lower than 
EDisCS. This is an indication that for the largest galaxies the 
majority of the morphological transformations took place a 
few billion years ago, while for most of the compact galaxies 
both the quenching of star formation and the final morpholog- 
ical type were reached at earlier epochs. 

It is clear that when comparing high- with low-z samples, 
it is of paramount importance to keep in mind that morpho- 
logically selecting galaxies at different epochs introduces an 
apparent, but spurious size evolution with redshift, which in- 



stead is a selection effect. Although it is impossible from the 
WINGS data to recover which galaxies were early-types at 
the EDisCS's epoch, our findings support the hypothesis that 
the main reason why the median size of WINGS early-type 
galaxies (dashed line in Figure |2|i is much more consistent 
with the median size of all EDisCS galaxies (dotted line, see 
also Table |2]l than with the size of only EDisCS early-types 
is that the largest late-type EDisCS's galaxies have gradually 
become earlier types by the WINGS epoch. 

5. DISCUSSION AND CONCLUSION 

We have found that 41% of EDisCS galaxies with M* ~ 
4x1O'"M0 are SDGs. Their properties are similar to WINGS 
SDGs, apart for a significantly different morphological mix: 
the prevalence of SOs in WINGS is not found in EDisCS. 

Such a result is not unexpected, given our previous findings: 
in VIO we have found that 17% (for the mass limits and radii 
adopted here) of WINGS clusters members at z ~ are SDGs. 
More than 50% of them have stellar ages older than 9Gyr, 
a clear indication that they were already old and compact at 
the EDisCS's epoch. The evolution of the SDG fraction in 
clusters with redshift is expected if SDGs are massive and old 
galaxies, formed in clusters seeds and preferentially found in 
today's massive clusters, while they are rarer in the field (see 
iTaylor et al.ll2009l) and therefore in the population of galaxies 
infalling into clusters at later and later times. 

We find that when galaxies of all morphological types are 
considered, the median size of cluster galaxies at z ^ 0.7 is 
only a factor 1.18 smaller than the local median. We conclude 
that from z ^ 0.7 to z ^ 0.04, there is at most a very modest 
evolution in galaxy sizes in clusters. 

Similarly to our VIO analysis of age selection effects, we 
have shown that comparing high-z morphologically selected 
samples with local ones can be misleading. In agreement 
with previous results regarding the morphological evolution 
in clusters, we have found that the largest EDisCS late-type 
galaxies are found to be large early-types in WINGS clusters, 
as it is apparent studying the morphologies above the tilted 
line in Figure|2]The BCGs, instead, have been found to evolve 
both in mass (a factor of ^ 2) and size (a factor of ~ 4), in 
agreement with other recent theoretical and observational re- 
sults. 

Our findings show that the progenitor bias (in age or mor- 
phology) plays an important role in the size-growth paradigm, 
and must be carefully taken into account when comparing lo- 
cal galaxy sizes with those of massive high-z galaxies. 
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